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I.  INTRODUCTION 

The  fact  that  acoustic  waves  could  give  optical  effects  in  nematic  liquid 

crystals  has  been  known  since  1936.  In  that  year  Fredericks  and  Zolin*  excited 

a nematic  liquid  crystal  by  tuning  forks  of  200  to  600  Hz  while  making  optical 
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observations  with  crossed  polarizers.  It  was  not  until  1969  that  Fergason 
suggested  that  a liquid  crystal  device  may  be  made  for  detecting  acoustic  in- 
tensity. Since  1970  there  have  been  approximately  forty  research  articles  pub- 
lished in  this  area.  However,  the  basic  mechanism  causing  the  effect  was  not 

3 

understood  until  the  last  three  years.  Differential  acoustic  absorption, 
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transverse  second-order  stress,  cybotactic  groups,  the  piezoelectric  effect 

and  anisotropy  of  acoustic  speed  have  all  been  proposed  to  be  the  mechanism 

for  the  effect.  In  the  early  work,  both  experimental  and  theoretical,  a 

threshold  of  acoustic  intensity  was  reported  to  be  required.  In  1976,  three 
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groups:  Sripaipan,  Hayes,  Fang;  Candau,  Peters,  Nagai;  and  Dion,  DeForest; 

reported  that  there  is  in  fact  no  threshold  for  the  effect.  The  former  two 

groups  attributed  the  effect  to  acoustic  streaming.  Further,  whereas  all  of 

the  previous  work  had  required  some  special  property  of  a liquid  crystal  for 

the  effect  acoustic  streaming  would  result  in  any  viscous  liquid.  It  is  simply 

that  acoustic  streaming  or  any  flow  in  a liquid  crystal  results  in  optical 

effects.  The  streaming  theory  has  been  extended  to  include  the  simultaneous 

effect  of  an  electric  field  by  Hayes*0  in  1978.  Again  the  streaming  model 

successfully  explained  the  experimental  results.  The  logical  result  of  the 

streaming  explanation  is  that  a cell  to  visualize  an  acoustic  field  could  be 

made  by  separating  the  liquid  crystal  into  distinct  regions  to  confine  the 
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flow.  Nagai  and  Iizuka  ’ have  made  such  a device.  However,  there  has  not 
been  a direct  observation  of  the  flows  induced  in  these  cells  nor  quantitative 
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theoretical  work  performed  giving  the  speeds  of  flow  or  the  flow  patterns  in 
the  cell.  It  is  the  goal  of  this  report  to  provide  direct  evidence  for  the 
flows.  In  Section  II  a theory  is  presented  which  predicts  the  magnitude  of 
the  flow  speed  and  the  form  of  the  flow  pattern  one  would  expect  in  a vis- 
cous fluid.  In  Section  III  results  are  presented  for  observations  of  these 


fluid  flows. 
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II.  THEORY  OF  FLUID  MOTION 
A.  Non-viscous  Limit 

We  will  develop  equations  describing  the  motion  induced  in  a disc  shaped 
liquid  bound  between  two  solids.  The  solid  surfaces  are  assumed  planar  with  an 
ultrasonic  wave  traveling  from  one  solid,  through  the  fluid  and  on  through  the 
second  solid.  We  will  take  the  cylindrical  boundary  of  the  disc  to  be  open  to 
a gas  of  negligible  density.  Such  a system  could  be  made  by  placing  a drop  of 
water  between  two  glass  microscope  slides.  We  will  further  take  the  distance 
between  the  planes,  the  height  of  the  disc,  to  be  smaller  than  the  wavelength 
of  the  ultrasonic  wave.  The  oscillatory  motion  of  the  wave  will  produce  a 
static  flow  in  the  fluid,  acoustic  streaming.  We  will  find  the  form  of  the 
resulting  flow  patterns  and  the  magnitude  of  the  fluid  speed. 

In  order  to  understand  the  approximations  which  must  be  made  to  obtain 
a tractable  solution  we  will  first  investigate  a simpler  problem  where  vis- 
cosity is  omitted  and  the  disc  is  assumed  of  infinite  height. 

Consider  a column  of  non-viscous  fluid  with  equilibrium  density  py 
supported  by  a piston.  See  Figure  1.  We  assume  the  piston-fluid  interface 
is  located  at  z = 0 and  take  the  radius  of  the  fluid  to  be  R.  The  result  of 
the  piston's  velocity,  Ae  *l0t,  is  a compression  wave  which  propagates  along 
the  cylindrical  axis  of  the  column,  the  positive  z axis.  Because  of  the 
boundary  at  r - R there  will  also  be  laterally  induced  waves  as  well.  The 
equation  describing  the  pressure,  P,  in  the  column  is 


V2  P - ^ 


1 

c2  at7 


a) 


where  C is  the  speed  of  the  wave  in  the  fluid.  At  r - R we  must  have  F • 0. 
The  solution  to  Eq.  1 satisfying  this  boundary  condition  is 
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P - E A J (k  r)  e1  kn  2 " ia,t 
n n o n 


where  J ( ) is  the  Bessel  function  of  order  zero. 


k2  * £L  _ k"2 

n C*  n 


aR  = kn  R - 2.405,  5.520,  8.654,  11.792, 


where 


J (a  ) *=  0 

o n 


n = 0,  1,  2,  ... 


The  partial  differential  equation  relating  the  pressure  to  v,  the 
fluid  velocity  is  given  by 


pofr  ‘ - 7P 


We  have  at  z = 0 the  further  boundary  condition: 


v = A e 
z 


Combining  Eq's.  2,  6,  and  7 we  have 


p a)  A = E A k"  J (k  r) 
o nnno  n 


Operating  on  both  sides  of  Eq.  8 with 


we  obtain 


2Pn  to  A 


n °n  Ki  (“n^ 

Mow  solving  Eq.  6 for  vf  and  taking  r = R we  find 


v = z M ei  kn  2 " iwt 

R n RkjJ  e 
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From  Eq.  3 and  4 we  see  the  value  of  R may  be  adjusted  with  respect  to 
A to  give  k^  ■ 0.  This  condition  gives  a resonance  and  may  occur  for  modes 
n ■ 0,  1,  2,  ...  . The  corresponding  values  of  R are: 


R - .383  A,  .879  A,  1.38  A,  1.88  A,  . . . 


(11) 


Figure  2 shows  the  amplitudes  of  vR/A  at  z * 0 as  a function  of  R/A  for  modes 
n ■ 0 and  1. 

From  Eq.  3 we  see  if  R < .383  A then  k^  will  be  imaginary  making  the  wave 
damp  out  in  the  z direction  as  seen  from  Eq.  2.  The  amount  of  damping  depends 
on  how  close  to  resonance  the  radius  is.  In  Figure  3 vR/A  is  graphed  as 
a function  of  Z/A.  We  see  for  the  mode  0 the  damping  is  less  the  closer  R/A 
is  to  the  resonant  value.  We  also  see  that  the  mode  n « 1 for  R/A  ■ .38  damps 
out  more  quickly  than  mode  n - 0 even  for  R/A  ns  low  as  .22,  about  40Z  lower 
than  the  resonant  condition  for  that  mode.  We  also  see  from  Figures  1 and  2 
that  although  the  disc  is  being  driven  by  a piston  oscillating  in  the  z direc- 
tion the  major  motion  can  be  radial. 

From  Eq.  6 we  can  also  find  the  radial  velocity  as  a function  of  r: 


£ 2iA  J)  (qnr/R)  ei(k»  z~0>t) 
n ' * Rkj;  J,  (on) 


(12) 


From  Eq.  12  we  see  at  r = 0 all  modus  are  zero.  For  the  n m 0 mode  the  maximum 
velocity  occurs  at  about  75*  of  the  maximum  radius.  A different  situation  would 
exist  for  an  infinite  strip  of  fluid  between  two  planes.  In  that  case  the 
Bessel  functions  would  be  replaced  by  sin's  and  cos's  and  the  maximum  speed 
would  occur  at  the  fluid  edge.  For  the  n *=  1 mode  of  oscillation  the  fluid 
for  r < . 7R  is  out  of  phase  with  the  outer  section.  Again  the  maximum  speed 
for  this  mode  does  not  occur  at  the  edge  but  at  about  r ■ .35R. 
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B.  Viscous  Model 

With  the  aid  of  the  results  we  have  obtained  from  Section  I1A  we  may 
now  proceed  with  a model  closer  to  the  physical  system  under  study.  Rather 
than  an  infinite  column  for  the  fluid  we  will  assume  a height  H with  the  origin 
of  Figure  1 moved  along  the  2 axis  so  the  upper  and  lower  boundaries  occur  at 
2 “ H/2  and  -H/2  respectively.  Since  we  are  now  including  viscosity  we  will 
start  with  the  Navier-Stokes  equation: 

Po  [fj  + (y  * * -c2  Vp  + nV2  v + (C+-|n  ) V (V  • v)  (13) 


Since  we  are  assuming  the  hydrodynamic  variables  may  be  expanded  about  their 
equilibrium  values  we  will  for  now  drop  the  higher  order  effects.  Wc  will 
take  for  a rectilinear  coordinate  system  in  component  form: 

p0  Ir1  = -c2  + ny2vi+  (?+In>a!l  07 (14) 

If  we  operate  on  Eq.  2 with  9/9x^  and  sum  over  i we  have 
9V  • "*  -v 

p0  —j~  - -c2  v2  p + n v2  v • v + a +j  n)  v2  (V  • v)  (15) 

The  hydrodynamic  variables  roust  also  obey  the  continuity  equation? 

1^  + V • (p  V ) = 0 (16) 

where  again  to  first  order  we  have 

!f+  PoV  = 0 (17) 

Combining  F.q's.  15  and  17  and  assuming  a time  dependence  eiwt  we  have 


where 


Aq  V2p  = p 


A 
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itt  + *D)  _ c2 
wp0  1 to7 


(18) 


(19) 
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The  solution  to  Eq.  18  is 

P ■ £ e ” + A2n  JQ  O^r)  0 J 


iwt 


(20) 


where  we  have  dropped  the  i from  in  front  of  the  kj^  factor  since  we  will  assume 
R < .183  A and  from  Section  1A  wc  know  for  that  case  the  z dependence  is  real. 
We  also  have  analogous  to  Kq.  3 

k2  - k"2  - | (2D 

Ao 

In  order  to  simplify  the  solution  of  the  Navier-Stokes  equation  we  define  two 
new  functions  according  to 


and 


Now  we  have 


3<h  ^ 

V ■■  v1  + V1 

r 3r  9z 

v«  ■ 1i  - r Sr  <rl,,) 


V • v - V7  4>  ; 


(22) 


(23) 


(24) 


so  from  Fq's.  17  and  18  we  find  <|>  may  be  related  to  p: 

* . - 

Po 


The  z component  of  Eq.  13  to  first  order  is 


Po  - -C7  *£  + n V7  v2  + 
Using  Kq's.  17,  18,  19,  and  22  we  obtain 


Ll  \ 3(V  • 

V^'V  -a*- 


v) 


( 2 *> ) 


(26) 


where 


V7  C - f, 

iu»p0 


* - | & (rt) 


(27) 


k7  - k’7  - -i^Pa 
n n n 


Taking 


(28) 


the  solution  of  Eq.  27  we  take  to  be 


5-s[A.„J»(knr>ek4Z  +*.n'MV)e"kS2] 


-kn  2 1 i0Jt 

e 


(29) 


From  Eq.  28  then 


»[a 


kAz 


<>  ■=  E I A J.  ( k r)  e n + A 

n I 3n  1 n 


^„n  Jo  (knr)  e knZ"|  eiWt/k 
‘‘“on  J n 


(30) 


Eq's.  20,  25  and  30  may  be  used  with  Eq's.  22  and  23  respectively  to  obtain 


+ Jin 

k 


= iMApltp  [ k\[z  -kj{z] 

'r  n < p0  LAm  + Azn e J 


* [A3n  ^ - \n  e_k"Z]}  J,  <knr)f 


(31) 


ia)t 


and 


v = l [ A eKz  _ A e k"z] 

z n ^ p0  L m 2n  J 

-[A./"2  + \/ki2]}Jo<V>el“t 


(32) 


Assuming  a no 
experiment  we 
damping  s hlch 
the  disc.  We 
approximation 
we  find 


A,  " 
in  2n 


■slip  condition  at  z = +H/2  we  must  have  v =0  there.  For  the 

— r 

have  performed  H/ X = .055.  From  Figure  3 we  see  the  amount  of 
has  occurred  as  the  wave  traverses  the  vertical  distance  H of 
will  therefore  take  as  an  acceptable  although  not  excellent 
that  v^  = Ae  at  z = +H/2.  For  these  boundary  conditions 


“2Po  kn  A 


iwAk"  ctnJ  (a  ) (2k'  - k2  H) 
n 11  i n n n 


(33) 


2 k2  HA 

L3n  = A-n  = a J.  (a  ) (2k'  - k2H)  ekl['H/2 
n * n n n 


(34) 


For  the  case  under  study  k'  >>  kzH  and  k"H  <<  1.  If  we  assume  we  are  near  the 
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mode  0 resonance  condition  from  Figures  2 and  3 we  see  we  may  drop  the  higher 
modes.  Under  these  approximations  we  have 


Akc 


r a J (a  ) 
o i o 


-2z  + idW  (ek°z  - e~koz)  1 J (k-r) 


iijjt 


l o 


(35) 


2 A J0(kor)e 
a J (a  ) 

o i o 


iO)t 


(36) 


From  Eq.  28  we  see  that  k^  is  complex.  For  what  actually  occurs  is  a 

cos  uit  rather  than  the  exponential.  For  the  radial  component: 


Akr 


v - r-“-7 — r s - 2z  + 

r a J ta  ) ) 

o i o'  { 
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Reko  (z-H/2) 


cos  [Imk  (z-H/2)  + wt] 
o 


-Reko(z  + H/2)  , , 1 

-e  cos  [Imk  (z  + H/2)  + tot ] j > x J (k  r) 

o Jl  l o 


By  setting  r = R in  the  real  form  of  v^.  we  may  obtain  a value  of  the  radial 
amplitude  of  oscillation  compared  with  that  of  the  piston.  For  our  sample 
we  have  H = .032  cm,  Re  k'  = Im  k*  = 1112/cra  and  at  resonance  R = .22  cm.  The 
resulting  ratio  of  radial  velocity  to  maximum  piston  velocity  (which  equals 
also  the  ratio  of  radial  distance  amplitude  to  vertical  piston  distance  ampli- 
tude) is  plotted  as  a function  of  z in  Figure  4.  The  resulting  maximum  ratio 
shows  for  our  case  the  radial  motion  has  an  amplitude  less  than  20%  that  of 
for  the  piston.  Therefore,  the  infinite  ratio  for  mode  0 shown  in  Figure  2 
turns  out  to  be  of  the  order  of  .1  when  viscous  effects  are  considered.  Also 
in  Figure  4 a diagram  is  shown  which  depicts  the  deformation  of  the  disc  cross 


u^m 


section.  As  the  planes  move  upward  the  fluid  bulges  out  at  the  bottom  and  in 


at  the  top.  As  the  planes  then  move  downward  the  opposite  occurs. 
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Let  us  now  return  to  Eq.  13  and  examine  the  term  we  omitted.  In  light  of 
the  solutions,  Eq's.  35  and  36,  we  see  that  each  term  of  Eq.  13  oscillates 
sinusoidally  in  time.  Although  in  Eq's.  35  and  36  we  have  used  the  complex 
exponential  form,  since  the  velocity  is  real  the  term  we  earlier  omitted  from 
Eq.  13,  PQ  v • Vv,  will  have  a factor  of  cos^  tot.  Since  cos^  <ot  = (1  + cos  2tot)/2 
if  the  time  average  of  the  omitted  term  is  taken  a non-zero  static  value  re- 
sults. Therefore,  the  time  average  of  the  other  terms  in  the  equation  must 
also  give  non-zero  static  values.  Since  these  other  terms  involve  only  the 
hydrodynamic  variables  to  the  first  power  we  conclude  a static  value  must  be 
added  to  the  values  we  obtained  in  Eq's.  35  and  36.  This  static  flow  is  called 
acoustic  streaming.  We  will  denote  the  static  addition  to  the  variables  with 
a subscript  2 signifying  the  second  order  solution:  v^2,  P2  > vz2'  The  time 
average  of  Eq . 13  to  second  order  is 

p <v  • Vvr>  = -cJ  Ji+nV2  v2r  + U+-|n)§^  (V  • v2)  (37) 

p <v  * Vvz>  = -c2  nV2v2z  + U+^n)f^  (V  • v2)  (38) 

where  the  brackets  denote  a time  average  and  merely  involve  replacing  the 
2 

cos  wt  factor  with  1/2.  If  we  operate  on  Eq's.  37  and  38  by  3/3z  and  3/3r, 
respectively,  subtracting  the  results  we  obtain 


p < — (v  • Vv^)  - -^z  (v  • Vvz)>  = n I- 


3z  '*  ’V  3r 

Again  invoking  the  approximations  used  to  obtain  Eq’s.  35  and  36  along  with 
kMl  >>  1 we  find  we  may  identify  V1*  with  3l,/3z'1  and  obtain  with  the  help  of 


(39) 


Eq's.  22  and  23 


p A*N^  Ji(knr)  < 3 J ] (knr) 

2na?  J (a)  ’ 3 (k  r) 

O 1 o o 


4 (R  k’)3 
e o 


3. 

c (cos 


a - sin  a)  + 


-b, 

c (cos 


h + sin  b)  + 


CjZ^  + CjZ^'b 


C3Z 


+ C. 


(40) 


n 


I •••  m a s~  /■  ■ 7- — ^ i H 
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where 


a = Re  k'  (z  - H/2) 
o 


b *=  Re  k'  (z  + H/2) 
o 


(41) 

(42) 


and  where  we  have  again  taken  the  real  part  assuming  Re  k^  = Im  k^  as  we  find 

to  be  approximately  true  for  our  systems.  The  contribution  from  <J>  can  be 

2 

shown  to  be  small  so  that  Eq's.  22  and  23  become 


v2r  = 9|/V9z 

v„  = - — v-  (r^  ) 
2z  r 9r  2 


(43) 


The  values  of  the  C^'s  may  be  found  using  the  boundary  conditions  that  both 
and  V2Z  equal  zero  at  z = + H/2.  Finally  we  obtain  the  second  order 
contributions  to  the  velocity: 


= _ P A2  k*  H J,(knr)  . 9 J i (k0r) 


2r 


q a2  Jj (a0) 


9(k  r)  1 2 (Re  k 

o ( o 


L 


ea  sin  a + e b sin  b 


3zJ 


2 (Re  k')3  H2  8 (Re  k')3 
o o 


(44) 


7z 


I JL  (T  J (k  r)  8 Ji(kor)\ 

r»  a 2 J2(c0)  r 9r  \r  Vo'  9 (kQr)  ) 


j - 4 (Re^1 ) 3 [ea(cos  a “ sin  a)  + e b(cos  b + sin  b)J 


3z 


8 (Re  k')3  2 (Re  k')3  H2 

o o 


(45) 


Eq's.  44  and  45  will  be  the  basis  of  a comparison  with  the  experimental  results 


we  have  obtained. 


III.  EXPERIMENT 


A barium  titanate  crystal  was  used  to  produce  the  ultrasonic  wave  (cw) 
in  the  sample  under  study.  The  wave  was  transmitted  from  the  crystal  through 
the  conducting  glass  electrode  immediately  above  the  crystal,  through  vacuum 
grease  which  was  used  to  mount  the  sample  holder,  through  the  bottom  glass 
slide, through  the  sample  fluid  disc,  and  on  through  the  upper  glass  slide. 
Various  fluids  were  used:  water,  a cholesteric  liquid  crystal,  a nematic  li- 
quid crystal  and  a transparent  vegetable  oil.  The  oil's  viscosity  is  closer 
to  the  nematic  whose  motion  we  will  ultimately  consider.  Therefore,  measure- 
ments reported  here  concern  the  velocity  induced  in  the  oil.  To  measure  the 
magnitude  of  the  flow  in  the  disc  the  oil  was  doped  with  a hydrocarbon  base 
ferrofluid  (0.1%  concentration).  Aggregates  of  magnetite  can  then 
be  seen  in  the  oil  by  means  of  a polarizing  microscope.  By  measuring  the 
speed  of  these  particles  the  speed  of  the  oil  may  be  determined.  The  parti- 
cles were  of  the  order  of  one  micron  in  diameter.  The  thickness  of  the  disc 
was  .032  cm  so  the  particle  size  was  small  compared  to  the  distance  over  which 
the  speed  of  the  fluid  changed  appreciably.  The  diamter  of  the  disc  was  from 
one  to  four  mm,  small  compared  to  the  12.7  mm  diameter  of  the  barium  titanate 
crystal.  The  microscope  was  calibrated  so  the  vertical  position  of  the  par- 
ticle whose  speed  was  being  measured  could  be  determined  when  the  microscope 
was  focused  on  that  particle.  A sealed  container  was  also  used  so  that  the 
ultrasonic  wave  would  first  travel  through  one  cm  of  water  before  reaching 
the  fluid  disc.  Using  a capillary  viscometer  the  ratio  of  viscosity  to  den- 
sity of  the  oil  was  deterntioed,  n/p  = 0.66  + 0.03  cm^/sec.  From  F.q.  28  we 
sec  can  now  be  determined  if  we  know  the  frequency  and  wave  speed.  We 

take  c “ 1.46  x 10^  cm/sec  and  f = .26  MHz.  We  find  the  real  part  of 

small  so  k'  = 1112  (1  + i)/cm. 
o 
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From  Eq.  44  we  see  that  the  radial  flow  goes  to  zero  at  r ■ 0 due  to  the 

J (k  r>  factor.  It  becomes  larger  for  larger  values  of  r and  then  goes  to  zero 

near  the  edge  r = R due  to  the  derivative  of  J. (k  r) . The  value  of  v~  on  the 

1 o i z 

other  hand  Is  larger  at  r = 0 and  decreases  as  r is  increased  and  finally  is 
larger  but  with  a different  sign  near  r = R.  A diagram  of  the  flows  our  equa- 
tions predict  is  shown  in  Figure  5.  The  upper  diagram,  I,  depicts  the  flow 
pattern  although  actually  H is  much  smaller  with  respect  to  R than  we  have 
drawn.  The  lower  pattern  is  not  what  we  predict  from  our  equations  but  the 
flow  patterns  that  would  pertain  if  the  models  proposed  by  other  authors  were 
applicable.  If  the  fluid  is  observed  at  r = R/2  as  one  focuses  from  the  top 
to  the  bottom  the  particles  can  be  seen  to  flow  toward  the  center  near  the 
top,  toward  the  edge  near  the  middle  and  toward  the  center  near  the  bottom,  in 
accordance  with  the  pattern  I which  our  equations  predict.  A more  quantitative 
comparison  is  presented  in  Figure  6.  The  graph  shows  the  fluid  speed  in  the 
radial  direction  (toward  the  center  is  taken  as  positive)  as  a function  of  z. 

For  the  data  shown  R/A  = .21  and  the  frequency  was  .26  MHz.  The  taller  (labeled 
I)  peaked  curve  in  Figure  6 is  the  best  fit  of  the  quantity  in  brackets  in 
Eq.  44  using  our  experimentally  determined  value  for  k'^.  Since  everything  in 

the  brackets  is  determined  a best  fit  was  made  of  the  data  to  a constant  times 

9 3 

the  bracket.  The  constant  so  determined  was  (13  + 1)  x 10  cm  sec/micron  and 
2 

X = 2.3.  The  smaller  peaked  curve  (labeled  II)  in  Figure  6 is  a best  fit 
where  both  k'^  and  the  constant  arc  to  be  found.  Their  values  are  respectively 
470  + 10/cm  and  2.0  + o.l  x 10  cnTyCec /micron  with  X = 0.3.  The  later  value 
of  k'  appears  unacceptable  in  light  of  the.  small  error  on  viscosity  and  the 
better  fit  we  take  to  be  fortuitous.  It  is  noteworthy  that  the  fit  is  so  good 
since  for  the  data  plotted  R/A  is  about  half  of  the  resonant  value  for  which 
Eq.  44  was  derived.  Notice  the  peak  near  the  top  of  the  disc,  the  one  on  the 


right  in  Figure  6,  is  slightly  less  than  the  one  on  the  left.  The  reason  for 
the  difference  has  to  do  with  our  assumption  that  v ■ Ae*a,t:  both  at  the  top 
and  bottom  of  the  disc.  Making  the  boundary  condition  that  = 0 at  the  top 
decreases  the  top  peak  to  about  20 7,  of  the  bottom  peak.  The  data  seems  to  be 
in  between  these  two  extremes  and  it  would  appear  we  have  taken  the  better 
approximation  in  allowing  the  upper  boundary  to  move.  Data  was  taken  for  R/X 
equal  to  .21,  .29,  .30,  .36,  .39,  and  .45.  All  of  the  radial  velocity  pro- 
files have  the  same  basic  pattern  as  those  data  shown  in  Figure  6. 

Frequencies  of  .26,  .77,  .68,  and  1.07  Mllz  were  used  both  for  direct  con- 
tact of  the  sample  holder  with  the  barium  titanate  crystal  and  with  the  crystal 
in  a cell  so  that  the  wave  was  propagated  through  about  one  cm  of  water  before 
reaching  the  sample.  All  of  the  flow  patterns  observed  in  the  microscope  using 
the  water  cell  were  radial.  For  the  direct  contact  case  for  certain  frequen- 
cies and  locations  on  the  crystal  it  was  possible  to  induce  patterns  of  either 
two  or  four  rather  circular  flow  patterns.  The  discs  became  slightly  distorted 
from  circular  to  ellipsoidal  in  shape  for  these  non-radial  patterns.  These 
patterns  are  assumed  due  to  non-uniform  oscillation  of  the  crystal.  For  these 
cases  the  pressure  in  Eq.  1 must  be  assumed  to  be  a function  of  angle.  The 

solutions  are  similar  to  the  well  studied  oscillating  drumhead. 

2 

The  radial  velocity  in  Eq.  44  is  proportional  to  A , the  square  of  the 
piston  velocity  amplitude.  But  this  amplitude  should  be  proportional  to  the 
voltage  applied  to  the  barium  titanate  crystal.  In  Figure  7 therefore  we  have 
plotted  the  radial  velocity  measured  at  a fixed  r and  z as  a function  of  vol- 
tage. The  solid  line  is  a fit  to  a voltage  square  curve. 
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Figure  1.  Diagram  for  calculations  of  the  infinite  non-viscous 

cylindrical  fluid.  The  piston  driving  the  oscillations 
is  in  the  xy  plane.  For  the  viscous  disc  of  height  H 
the  origin  is  moved  along  the  z axis  a distance  H/2. 


Figure  2.  Graph  for  the  ratio  of  radial  speed  amplitude  to  piston 
speed  amplitude  versus  the  ratio  of  cylindrical  radius 


to  acoustic  wavelength 


Figure  3.  Graph  of  velocity  ratio  versus  z/\.  The  velocity  ratio 
is  the  ratio  of  radial  speed  amplitude  to  piston  speed 
amplitude. 
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Figure  A.  Velocity  ratio  (radial  to  piston's)  versus  distance 


along  the  z axis.  The  model  assumes  no-slip  at  the 


boundaries  and  incorporates  viscosity. 

►. 


Figure  5.  Flow  patterns  in  the  disc  due  to 


The  upper  diagram,  I,  is  the  pattern  predicted  from 


our  second  order  solutions.  The  lower  pattern,  II 


is  that  assumed  in  other  theories 


Figure  6.  Velocity  of  fluid  due  to  acoustic  streaming  versus 
distance  through  the  disc  at  r = R/2.  Curve  I is 
a best  fit  to  the  data  using  our  measured  value  of  k' 
Curve  II  is  a best  fit  allowing  k'Q  to  be  determined 
by  the  fitting  process  itself. 
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Figure  7.  Velocity  of  the  fluid  due  to  acoustic  streaming  versus 
voltage  applied  to  the  transducer.  The  solid  line  is 
a best  fit  to  a voltage  square  dependence. 
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